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A model for classical impurities moving in Bose-Einstein Condensate (BEC) is proposed in the 
framework of quantum field theory and is solved within the Bogoliubov approximation at zero 
temperature. Several formulae are obtained for physical quantities such as the occupation number, 
the dissipation, the free energy, and the depletion. To illustrate this model two examples are studied; 
an impurity moving with i) a constant velocity and ii) a constant acceleration. Landau's criterion 
for energy dissipation accompanying the motion of impurities due to creation of quasi-particles in 
a superfluid, which was proposed on the basis of purely phenomenological argument, is obtained in 
the framework of an exact quantum mechanical model. 



The recent experimental discovery of BEC in alkali va- 
pors [1] was facilitated by magnetic trapping and the 
development of novel cooling techniques. This discovery 
and subsequent experiments motivated extensive theoret- 
ical activity and studies of BEC with additional trapping 
potentials [2]. 

Impurities are present in all experimental situations 
and should be accounted for. Here we discuss dissipation 
associated with motion of impurities in BEC. A semi- 
nal work on the problem of motion of impurities in a 
superfluid is due to Landau [3, 4]. He proposed a sim- 
ple criterion for energy dissipation based on purely phe- 
nomenological considerations demonstrating that an im- 
purity moving with a velocity greater than some critical 
value generally produces excitations in a fluid. Super- 
fluidity is explained then by the absence of excitations 
for motions with velocities smaller than the critical Lan- 
dau velocity. We sketch here shortly Landau's original 
argument for the sake of later comparison with the re- 
sults of our model computation of the effects of motion 
of impurities in the weakly degenerate Bose-Einstein gas. 

Suppose the excitation in a fluid, whose atoms have 
mass M, has an energy E and a momentum p in a ref- 
erence frame where the fluid is at rest. Then the en- 
ergy of the excitation E 1 in the coordinate system where 
the fluid is moving with velocity v with respect to the 
container, is given, in classical mechanics, by a Galilean 
boost: E' = E + p ■ v + Mv 2 /2. Dissipation due to the 
motion of the excitation can happen if the energy in the 
boosted frame is negative, E + p ■ v < 0. In the optimal 
case of p anti-parallel to v this implies v > E/p. There- 
fore, for the phonon-like excitation E — pc this condition 
gives v > c, with c the speed of sound in the fluid. While 
this condition is widely accepted, it seems that it was not 
derived yet from an underlying microscopic theory. 

Static "quenched" randomly distributed non-dynamic 
impurities were studied in [5, 6]. The problem of impu- 
rities moving in condensates has not been studied suffi- 
ciently in its full generality until relatively recent partial 
results [7-10]. Dynamical quantum impurities were also 
studied recently [11, 12]. 

The purpose of this Letter is to derive Landau's cri- 



terion from a microscopic quantum many body theory, 
and to study more general effects of impurities in BEC 
in a consistent and unified approach. To this end we 
propose a model Hamiltonian with classical impurities 
moving along given trajectories in BEC and diagonalize 
it using Bogoliubov's method [13]. Landau's criterion 
appears naturally. We also compute the effects of an ac- 
celerated impurity that can hopefully be experimentally 
tested. 

The Bogoliubov Hamiltonian for interacting bosons of 
mass M is 

Ho = J d 3 x [ft(x) (-^Pj fa) 

+ 9 -^{x)fa{x)fax)fax)}. (1) 

The coupling constant g is related to the s-wave scatter- 
ing length a s by g = Ana s h 2 /M. We suggest adding the 
interaction Hamiltonian; 

H- mt (t) = K [ d 3 x S(x - b(t))ft(x)i>(x), (2) 

i ■* 

to account for the local interaction between bosons and 
impurities moving along given trajectories Ci(t). The 
coupling constants Xi between i-th impurity and bosons 
are related to the corresponding s-wave scattering lengths 
b l s and reduced masses m, by Xi — 27tb\h 2 /rrii. Expand- 
ing the fields %j}>(x) and ip(x) in terms of creation and 
annihilation operators in the plane wave basis in a box 
of volume V — L 3 , we obtain, 

+ E#E^ (fe - fc ' Kl(t) 4^ (3) 

i k,k' 

where Ck — h 2 k 2 /2M is the free particle kinetic en- 
ergy. In the thermodynamic limit (N, V — > oo with 
number density n = N/V finite), limv_ >00 (l/y) J2k = 
(2tt)- 3 / d 3 k, lkm/^oo ^4 = (27r) 3 / 2 at(fc), and 



2 



liniy^oo Wd k = (27r) 3 / 2 a(fe). Then fit(fc) and fi(fc) sat- 
isfy the usual commutation relations: [fi(fc), o^(fe')] = 
5(fc - fc'). [fi(fc), fi(fc')] = [a f (fc), a f (fc')] = 0. 

In the Bogoliubov original treatment [13] of the above 
system the total number of particles is not conserved. In 
the following we adopt the standard approach [14-16]. 

We use the identity fijfio = N — ^2 o fc &fc to eliminate 
the zero mode number operator (a prime indicates the 
omission of the zero mode). In a large TV approximation 
we keep terms up to second order in at and fife. The 
Hamiltonian then commutes with N and choosing an ap- 
propriate basis replaces TV by its eigenvalue TV. and do 



can be replaced by \/iVo = y N — d\dk — VTV [17]. 

In order to diagonalize the Hamiltonian (3) we first 
implement the Bogoliubov transformation in the same 
approximation. The required canonical transformation 
Ub = exp(iGb(&)) is generated by 



G6(0) = ^ , M fe fi t _fe + ii.c. 



(4) 



with 



k = tanh 1 ( ^— ) , hcj k = y/e k (ek + 2gn). 

\nui k + e k + gnj 

(5) 

Bogoliubov's excitation (the "bogolon") is created and 
annihilated by 

b{ = e^Wfit e -«V«) = fit cosh0k + & _ k S mh6 k , (6) 

b k = e i6b ^a k e' ldb ^ = a k cosh0 fc + al k smh0 k . (7) 

The original vacuum state |0) is transformed to the new 
bogolon vacuum: 



\0(9)} = e l6b{9) \0) 
r' 1 



(8) 



n'e^ exp (^ tanMfe ) |o> - 

The Hamiltonian after this Bogoliubov transformation is 
H B (t) = E + J2'^kb{b k + Uk(t)b{ + h.c), (9) 

k k 

where the ground state energy Eq with no impurities is 
E a = gn 2 V/2 + J2'( h ^k -e k - gn)/2 = (gn 2 V/2)[l + 
(128/15) v /nof77r], and 



h{t) 



ne k 
Vhuj k 



e~ ik <^. 



(10) 



The Hamiltonian (9) describes a system of decoupled 
forced harmonic oscillators and, therefore, can be solved 
exactly [18, 19]. Diagonalization of (9) is easily accom- 
plished by way of time dependent unitary transformation 
Uc = exp (iG c (t)) generated by 



(11) 



This transformation shifts b\, and b k - 

e *<5e(*)gt e -iG c (t) =gt 

e iG c (% k e -iG c (t) = ft fc + fc ( t ). 



(12) 
(13) 



In order for G c (t) to diagonalize the transformed Hamil- 
tonian in the Schrodinger picture, 

H c (t) = e i6 ^H B {t) e- ld ^) _ lh e^d t e~ iC} ^\ 

(14) 

4>k(t) should satisfy a simple differential equation: 



ihd t 4>k{t) = huj k (j) k (t) + fk(t). 



(15) 



We solve this differential equation using the initial con- 
dition (f>k(to) — 0. It is assumed that the interaction is 
switched on at t = to- 



ne k 



-i(Aj k t 



h V Vhuj k 



f dt' 

i J t 



; iw fe t'— ik-Ci(t') 



(16) 



After the unitary transformation is performed we find 
that the Hamiltonian (14) is diagonalized as 

H c {t)=Eo{t)+Y!^kb{bk- (17) 



The energy spectrum hui k remains the same as bogolons, 
the gapless excitations characterized by u k ~ kc for a 
small k, where c = y 'gn/M is the speed of sound. In the 
Hamiltonian (17) we have defined a new ground state 
energy of the system by 



E (t) =Eo + lY! VhWMt) + c - c -), 



(18) 



with the second term representing the effects of impuri- 
ties. The bogolon vacuum is shifted corresponding to the 
unitary transformation !7c as 



|6(t)> 



iG s (t) 



10(f)) 



(19) 



ex P [-£'( hMt)\ 2 + Mt)b{W(e)). 



The physical interpretation of the vacuum (19) is quite 
simple. The motion of impurities creates time dependent 
coherent states. The new vacuum describes the deformed 
condensate and contains infinitely many bogolons [20]. 

We next proceed to evaluate the occupation number 
ilk from which various physical quantities such as the 
dissipated energy, the depletion of the condensate, etc. 
can be derived. 

The expectation value of the number operator h k — 
b k bk with respect to the shifted original vacuum 
exp(iG c (i))|0) is 

n k = <0|e-«*« b{b k e«*«|0> = fe*±g+ig +ft fc (t), 



2huj k 



(20) 
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where fife (t) = |0fc(t)| 2 is the contribution to n k due to 
the motion of impurities. From Eq. (16) we find 



ne k 



hk ^ = w^- lE A ^)l 2 < ( 21 ) 



where Ii(t) is defined by 

It(t) = I dt 1 e^t'-ifc-WO. 



(22) 



The energy £ k (t) dissipated due to the motion of impu- 
rities in a given mode k is calculated as the expectation 
value of the occupation number with respect to the new 
vacuum (19) multiplied by the exciation energy %u k , 



£ k (t) = {Q{t)\hu k b{b k \Q(t)). 



(23) 



The total dissipated energy £ t ot(£) is the sum over au 
modes fc: 



x — ^ f Ti x — x — ^ 

^tot(i) = 2^ huj k n k (t) = JT^2-^ 6k IZ^ A *^(*)I 



The depletion of the condensate is defined by 



(24) 



rf(*) = ^E'w*)i4w)>- (25) 

k 

The straightforward computation gives the result 



k 

-jnRe(^(t)^(i))]. (26) 

In the thermodynamic limit the second term is of the 
order of 1/N and will not contribute to the depletion. 

To illustrate our model consider the simplest case of 
one classical impurity moving with a constant velocity 
v = vz in the z-direction. For concreteness assume 
that the particle starts moving from one end of the box 
[z = 0) at t = and arrives at the other end of the 
box (z = L = y 1 / 3 ) at t = tf = L/v. Our main ob- 
jective is to calculate the occupation number fife (t) in 
the thermodynamic limit of large t. Therefore, we need 
to compute Jqo = li m t/->oo I(tf). We then face a formal 
mathematical difficulty because |ioo| 2 is the square of the 
distribution: 



lim 

7J->0+ 



dte iwt -^ =V(-)+it5(w), (27) 



where V is the principal value distribution, 5 is the Dirac 
distribution, and w — uj k — k ■ v. The standard way to 
deal with this difficulty is as follows, 



\I(t)\ 2 = 2 



1 — cos(wt) 



~ 27rf S(w), 



(28) 



where t is large. Thus |/(£)| 2 is proportional to t with a 
coefficient that becomes a delta function in the limit t — > 
oo. Lastly, we remark that the replacement \I(t)\ 2 /t — > 
2ir 5(w) is valid only asymptotically and the actual peak 
always has a finite width. 

We find that the rate of energy dissipation in the ther- 
modynamic limit is 



7 ■ 



d£tot (t) 
dt 



(2 



^ J d 3 ke k S(u k -k-v). (29) 



Thus 7 = unless there are such values of k for which 
u> k = k ■ v. To proceed further, we rewrite (29) in 
spherical coordinates and introduce new variables x = 
hk/(2Mc), y — f3cos9, where (3 — v/c. The result of 
evaluation of the integral (29) is 



7 = 



4nA 2 M 3 c 
irh 4 [3 



f dxx 3 f dy 5(Vx 2 + l-y). (30) 

Jo J-B 



Clearly it follows from Eq. (30) that f3 has to exceed 1 
in order to have energy dissipation and hence, Landau's 
criterion is derived. The rate of energy dissipation is: 



7 



n\ 2 M 3 c 3 



(/3 2 -l) 2 6(/3-l), 



(31) 



where Q(x) the step function. Like Cherenkov radia- 
tion this radiation is emitted inside a cone with opening 
angle 8 C — cos" 1 (l//3), and its momentum cuts off at 
hk c = 2Mc^ ' P 2 — 1 keeping the total dissipated energy 
finite. Our result agrees with that obtained previously in 
another approaches based on the Gross-Pitaevskii equa- 
tion [8-10]. 

As another example consider a single classical impu- 
rity moving with a constant acceleration a in the z- 
direction. We adopt the relativistic form of the tra- 
jectory C(t) = (c 2 /a)sj\ + (at/c) 2 z (which reduces to 
af 2 /2+constant in the non-relativistic limit) so that the 
velocity of the particle is always bounded by the speed 
of sound \C{t)\ < c. To obtain analytic results we extend 
the range of integration in (22) to (— oo, +oo), taking into 
account the complete trajectory. This corresponds to an 
idealized thermodynamic limit. In reality we always have 
to use finite times, volumes, and trajectories conforming 
to the experimental situations. 

The integral (22) is now evaluated, after the change of 
variables from t to r, t = (c/a) sinhr, 



a J- { 



dr e m sinh < T -x) coshr 



2c 

— (ir5(£l) coshx + iKi(£!) sinhx), 



(32) 



(33) 



where fl — (cja)^Juj\ — (ck ■ £) 2 , x = tanh _1 (cfc ■ z/u> k ), 
and Ki(x) is a modified Bessel function of the second 
kind. Since Q = if and only if k = 0, the first term in 



4 



Eq. (33) does not contribute to the occupation number 
after multiplying by e^/uk in Eq. (21) and using x 5{x) — 
0. Hence we obtain 



4c 2 n\ 2 e k . 2 

n k = 77Z* n — ( K i("J sinhx) ■ 



Vn 3 a 2 iv k 



(34) 



In contrast to the constant velocity case there is an ex- 
ponential tail at high momenta (which follows from the 
asymptotic behavior of Ki(x) for x — > oo), 



167rnA 2 M 3 c 3 cos 2 ( tick 2 
n k - , , exp - 



Vh 5 ak 4 



Ma 



(35) 



where 9 is an angle between the directions of the momen- 
tum Tik and the trajectory of the impurity. The formula 
(35) displays the exponential dependence on the mass M 
of a heavy boson. The particle count by a detector can 
be used to measure M. The total dissipated energy is 
evaluated in the weak acceleration limit by using again 
the asymptotic form of Ki(x), 



£tot — 



nA 2 M 2 c / ha 



ttMc 3 ' 



(36) 



Thus we find that although the velocity of the impurity 
never exceeds the speed of sound, there is a finite amount 
of energy dissipated due to the motion of impurity. This 
result may appear couterintuitive. However this is not 
an artifact due to the presence of an infinite time inter- 
val in our calculation. It simply means that Landau's 
phenomenological argument based on the Galilean co- 
variance does not apply to the constant acceleration case. 

As we have seen, our simple model Hamiltonian (3) 
yields Landau's criterion in the case where the impurity 
moves with a constant velocity. In contrast to the original 
argument, however, we find finite energy dissipation even 
for slightly sub-critical velocity v. This can be seen from 
Eqs. (28-31). 

We have also obtained the dissipated energy for a 
constantly accelerated impurity in the weak acceleration 
limit. In order to estimate experimental dissipation for 
this case and the previous case, we should evaluate ex- 
pressions like (24) in a way that corresponds to the actual 
numerical parameters. 

Our model is applicable to more general cases such 
as several impurities moving with different velocities or 
in a completely random way. Our model should be also 
extended to finite temperatures in order to deal with real 
experiments. 

Finally, we point out that the exact expression for the 
occupation number (34) shows that the vacuum conden- 
sate does not mimic a simple thermal bath usually as- 
sociated with the accelerated motion of impurities. The 



results of an extended analysis of the latter problem will 
be presented in the forthcoming publication [21]. 
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